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TECHNICAL NOTE 3525 

VORGZEX INTERFERENCE ON SLENDER AIRRLANES 
By Alvin H. Sacks 


SUMMARY 


Formulas are developed for the forces and moments (except drag) due 
to vortex interference on slender •^-ring-hody-tall combinations of general 
cross section in terms of the positions and strengths of the shed vortices. 
The analysis is applicable to steady motion and to motions which can be 
considered to be made up of a succession of steady states (i.e., quasi- 
steady motions). In order to illustrate the application of the analysis, 
the interference lift of a plane wing-body-tall combination in steady 
straight flight is determined by utilizing vortex positions obtained by 
numerical methods. 

It is found that the impulse of each shed vortex and its image 
vortex in a transformed circle plane enters into all the interference 
forces and moments on the airplane. A simple theorem is given for the 
interference forces in steady straight flight idiich are found to depend 
on this impulse evaluated only at the -I'/ing trailing edge and at the base 
of the configuration. 


INTRODUCTION 


It has been recognized for some time that interference among the 
varioiis airplane components can be very important in determining the aero- 
dynamic characteristics of slender configurations. The use of slender- 
body theory for treating entire wing-body combinations is not new (see, 
e.g., refs. 1, 2, and 3)7 for cases falling into this category con- 
siderable siniplificatlon has resulted. The calculation of interference 
effects due to wing wakes (wing-tail or wing-afterbody interference), on 
the other hand, is not so clear-cut. Unlike the iring-body problem, this 
calculation becomes more difficult as the effects become more important, 
since the rolling-up and displacement of the wing vortex sheets cannot 
be ignored for long slender configurations. 

Various authors have treated certain specific cases of interference 
due to VTing wakes by assuming simplified irake shapes (e.g., refs. 4 to 7). 
The present paper, however, is concerned with developing formulas for 
ceilculating the forces and moments due to wing wakes for slender wing- 
body-tail combinations; the formulas developed are in terms of the 
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strengths and positions of the shed vortices and the mapping fimctions 
of the airplane cross sections . Since the purpose of this report is the 
presentation of the derived formulas rather than the presentation of cal- 
culated results^ only the simplest example will he calculated here. 
Specifically, the interference lift of a plane wing-hody-tail combination 
iTlll he determined hy making use of a numerical analysis for the vortex 
strengths and positions at the tail. 


ANALYSIS 


It is generally recognized that in the vicinity of a slender hody 
the ^ra-ve equation in three space dimensions can he approximated hy the 
two-dimensional Laplace equation in planes normal to the hody axis, 
provided the frequency of the motion is small con 5 >ared with the flight 
velocity divided hy the length of the hody. Hence, the classical methods 
of hydrodynamics can he applied to quasi-steady motions of slender bodies 
in a compressible fluid. The analysis of reference 8 was the counterpart 
of a method due to H. Blasius for obtaining the forces and moments on a 
two-dimensional hody in an incon 5 )ressihle stream, and the hody cross 
section was mapped onto a circle in order to express the slender-hody 
stability derivatives in terms of the mapping functions of the cross- 
sectional shapes. The pressures on the hody were expressed in terms of 
the two-dimensional complex potential and the total forces and moments 
were calculated hy integrating those pressures round the hody cross 
section and over the hody length. The complex potential was expressed 
in a Laurent series of the form 


00 

F = B in ^ +2^ — + D 
n=i ^ 

and the integrations were carried out hy the method of residues . In the 
present analysis, it ^ ri 1 1 he shown that hy a simple extension of the 
resulting formulas of reference 8 for the forces and moments the effects 
of ■t'dng wakes can he included so that 't/lng-hody-tall combinations may he 
treated in the same fashion. 

The extension to he made icLll consist of admitting any number of free 
vortices external to the hody cross section along with their corresponding 
Images inside the hody. The boundary conditions are unchanged, then, and 
the only fundamental change lies in the expression for the complex poten- 
tial which enters into some of the contour integrals to he taken romd 
the hody cross sections for the determination of the forces and moments. 

It ^fill first he observed that all the vortices in the flow field 
external to the hody will he free vortices, representing wing wakes, and 
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ttierefore can themselves sustain no forces. Hence the integrals for the 
forces, ■vThen evaluated round each external vortex, must vanish. In vie'w 
of this fact, the resulting expressions for the forces and moments will 
he the same whether the integrals are evaluated at the body surface or 
on any other contour enclosing the body cross section. The integrEils 
involving the potential can therefore be evaluated at large distances 
from the body so that, just as in reference 8 , the con5)lex potential will 
he expressed as a Laurent series valid at large distances, the con- 
tour integrals will again he evaluated by the method of residues . 


The form of the complex potential will he exactly the same as eq^ua- 
tion (1) since the same types of singularities are introduced by the 
addition of external vortices and their images. Therefore, any contom* 
integrals which were identically zero before the addition of the vortices 
■vd.ll remain zero. Fin-thermore , it was found in reference 8 that all the 
forces and moments (except drag) are linear in the potential. Hence it 
is possible to calculate the forces and moments due to vortex interference 
alone (by using only the additional potential due to the vortices and 
their images) and then to add these directly to those calculated idthout 
interference. The present analysis vdll therefore he concerned -with the 
determination of the additional potential due to the shed vortices and 
their images and ^dth the associated interference forces and moments. 

If one considers a slender w ing - 
hody-'fcail combination as shoim in 
sketch (a), it is noted that the 
net circtilation in planes perpen- 
dicular to the X axis must be 
zero for all values of x. Hence 
if the airplane cross section in the 
physical plane ^ = y + iz is trans- 
formed to a circle of radius r© 

(center at the origin) in the 0 
plane by the transformation 



^ = f(u) = 0 +y 

n=o (2) 


then the additional complex poten- 
tial F’ due to the shed vortices 
and their images can be expressed 
in the plane of the circle simply 
as (see sketch (b)) 



( plane <r plane 


Sketch (b) 
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F’(^) = Fi»((t) = - - Ok) + -^^rk2n(a - 

k=i k=i 


(3) 


where m is the number of extemeO. (free) vortices and may he a function 
of X. All other symbols of the present analysis are defined in Appendix 
A. Now, since the transformation of equation (2) leaves the flow field 
at infinity unchanged, the residue of the complax potentieil is the same 
in either plane. Thus, expanding equation ( 3 ) for large 0 , one finds 



so that the residue Ax’ is seen to he 


Ax’ 


m 



k=i 


m 



k=i 


(5) 


where represents the distance between the kth free vortex and its 

imsLge in the complex 0 plane. It is interesting to note that the 
q\iantity fkok^ is proportional to the impulse of the vortex pair made 

up of the kth external vortex and its image in the 0 plane. 


According to reference 8 the total forces and moments on a slender 
configuration whose trailing edges all lie in the base plane x = Z are 
given by 


S(R + 2ipCc) + Uo ^ (S^,,) 


x=Z 


Y - iL = atpUoAi^^^ + pUo 

r ^ 81 ? r^— r ^ — — 

p / S ^ dx + 2jtipp / .Axdx + ipp / S(R + ip^(,)dx 


I 


8Ax 


+ 2itp / — dx + 


at 


( 6 ) 
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N - 121 = - 2 rtpUo / (x - Cl) dx - 


pU, 


/ s 


S(R + 2 ip^c) + Uo ^ (S^c) 


dx - 


I I _ 

2rtp J (x - Cl) dx - p y' (x - Ci)S ^ dx 


2jtipp y (x - ci)Aidx - ipp J' (x -Ci)(R + ip^£,)S dx 


( 7 ) 


L 2 

L« = I pUqRjT F d(?r) + I p^y dx ~J Fd(^I) - atpi y midx - 
x=Z o o 


ppRy S(RCc - 2R^c)<3^ - PUoI y (sCc)dx (8) 

o o 

where R denotes the real peirt, -I denotes the imaginary part, and p 
is the rate of roll about the body axis. Now, in order to write the 
corresponding expressions for the airplane of sketch (a), it will be 
convenient to divide the a±rplane length into the two segments shown. It 
should be mentioned that the present analysis will be restricted to I'/ings 
I'Tith traili n g edges normal to the x axis. However, if the wing trailing 
edges do not all lie in the same plane x = constant, further division 
trf.ll of course be necessary, but the procedure trfll be the same. Equa- 
tions ( 6 ), (7), and ( 8 ) can be applied directly to each segment of the 
airplane of sketch (a) to give the total forces and moments including 
vortex interference if we replace F by (F + F*) and Ax by (Ax + Ax*)^ 
where F* and Ax’ are the additional con^ilex potential due to the shed 
vortices (and their images) and the residue of that additional complex 
potential. Thus, making use of equation ( 5 ) for. the residue Ax’, and 
stipulating that the interference forces act only on portions of the 
airplane behind the T'rfng trailing edge, we find that the additional forces 
and moments due to vortex interference (and denoted by the subscript l) 
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wliere refers to a station immediately lehind tlie wing trailing edge 

Inasmuch as F’ has been expressed in the transformed circle plane 
(eq. (3)) it ^/ill, in general, he convenient to carry out the contour 
integrations of equation (U) in the a plane, particularly since the 
integrals become analytic in that plane due to the relation Ob' = ro^ 
on the circle hoimdary. 


For the special 
and (11) reduce to 


case of steady straight flight, equations ( 9 )^ (lO), 




Yj " ilq — “IpUo 


x=l 


( 12 ) 
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Nj - IMj = IpUo 



TE 




ix 


( 13 ) 


Li» = I pUoH jTp*d(^0 - I pUo» J 

X=2 X=XrjTg^. 


I ^ 

pR f S^rtot^dx (14) 
TE k=i 


m 


and it is interesting to note that the quantity 




enters the 


calculations in much the same manner as the apparent mass enters into the 
stability derivatives of wing -body combinations (see ref. 8). In fact, 
equation (12) can he stated in the form of a simple theorem for the 
determination of the Interference forces. 


1‘JHEOREM : The lateral force Yj + iLj due to each vortex of 

strength F shed from a fon^ard 'vring of a slender wing-hody- 
tail combination in steady straight flight is equal to the 
change, from -tiong trailing edge to base of the airplane, of 
the q-uantity ipUoTffi. where Op is the (complex) distance 
bet^reen the vortex and its image in the plane in which the 
body cross section is mapped onto a circle ■^diile leaving the 
flow field at infinity unchanged. 


An important point to notice in regard to equation (l2) n.nfl the 


above theorem is that the quantity 


m 

pUoB^ 


^Tk-r evaluated at the iring 


k=i 


trailing edge is equal to the lift of the airplane segment lying ahead of 
the -vring trailing edge. Hence the interference lift in steady straight 




flight is equal to the difference between the quantity pUoB.1 

\k=i /x=7 

and the lift of the airplane segment ahead cf the iirlng trailing edge. 
Thus, if we TOite for the total lift of the airplane of sketch (a) 


L — La + Lt) + Lj 

I’There Lg^ and L|^ are the lifts of the Isolated segments of the airplane 
lying ahead and behind the wing trailing edge (i.e., trLth no Interference), 
it is clear that 
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m \ 


L - La + Ljj + pTJosf 


k=i / 

^ '^x=Z 


- La. — L"b + pUqH, 



The total lift of the airplane is therefore equal to the lift of the 
isolated airplane segment behind the t-rLng trailing edge plus the quantity 

pTIqB/ ^ lift Ljj can often he calculated hy ordinary 


slender-hody theory and for many interesting configurations this part of 
the lift is already kno™. The above statements are independent of the 
form in I'/hich the vortex sheet leaves the wing, hut this will of course 
determine the positions and strengths of the vortices at the hase of 
the conf iginration . 


It can he seen from equation (12) and the associated theorem that in 
steady straight flight the interference side force and interference lift 
depend only on the shapes of the cross sections at the wing trailing edge 
and at the hase of the configuration and on the strengths and positions 
of the shed vortices at these two stations. This statement closely 
parallels the corresponding statement for slender -vring-hody combinations 
that the side force and lift depend only on the shape of the hase cross 
section (see ref. 1) and on the angle it makes 'v/lth the flight direction. 
The forces and moments including interference effects will, howver, 
he nonlinear with respect to the angles of attack and sideslip. 


It is evident from equations (9)j (lO) and (ll) that the calculation 
of the interference forces and moments requires, in all cases, a knowledge 
of the vortex strengths and positions as well as the mapping fimctions of 
the airplane cross sections. Q!he calcvilatlon of the forces and moments 
without vortex Interference requires only the latter. 


APPLICATION TO A PLANE NTNO-BODY-TAIL COMBIMTION 


z 



In order to illustrate 
the application of the fore- 
going analysis, a simple 
example iriLll he treated briefly. 
In this example the interfer- 
ence lift will he determined 
hy means of equation (12) . 
Consider a plane wlng-hody-tail 
combination in steady straight 
fll^t at an angle of attack a 
as shoT'/n in sketch ( c ) . The 
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cross section is taken as a circle with symmetrical flat-plate ■wrings 
mounted along a diameter. The complex variable is then given by 


cr 




and ^o = ^ where a 

is the body radius and s is 
the local semispan of the ivlng 
or tail. It vriLll he assimed 
that the trailing vortex sheet 
is fully rolled up into two 
voartices somewhere ahead of the 
tail trailing edge which lies 
in the plane x = I . Thus , 
denoting the station x = Z 
by the subscript 1 , we have 
at X = Z (see sketch (d)) 


z 



(X~l) (x-l) 


Sketch (d) 



Fin:thermore , since here by symmetry 


P2 = = -r 

and 

ta = -h 

it follows that 
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Now, from equation ( 12 ), 


Li = pUoR 



- PUqR 



and the second term haa already been identified as the lift of the air- 
plane segment ahead of the wing trailing edge. For the present problem, 
if the vortex sheet leaves the i^ing trailing edge as a flat sheet from 
each wing panel, this lift is known to be (see, e.g., ref. l) 



where aQ and Sq are the body radius and wing semispan at x = x^. 
Therefore, the interference lift for the plane -iring-body-tail combination 
of sketch (c) is 


- (3. . - 


jtpUo^aso^ 



and it is clear that F and (the strength and position of rolled-up 
vortex 1 at x = I) must be specified if the interference lift is to be 
calculated . 


Before any calcifLations are performed, it is interesting to note 
that if the real part of the sq\iare root in equation (15) vanishes, then 
the interference lift is equal and opposite to the lift of the airplane 
segment ahead of the wing trailing edge. This condition is satisfied if 



is real and less than 



and it can easily be shown 


that this requires that real and less than Sx. In other words, 

if the roUed-up vortices intersect the tail anywhere along its trailing 
edge, then the interference lift just cancels the lift forv/’erd of the 
vring tra iling edge. It can therefore be concluded that, for a I'ring-body- 
tail combination having cross sections of the type shown in sketch (d), 
if the location and span of the tail are such that the rolled-up vortices 
pass through the tail trailing edge, the total lift of the combination is 
unchanged by removing the portion of the airplane forward of the wing 
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trailing edge. This result has required no assumptions regarding the 
influence of the tail on the vortex paths hut has required the assumption 
of inviscid vortices (no cores). While the above result wotild be modified 
by the presence of viscous cores, equation (15) would still apply at angles 
of attack at which the vortex cores do not touch the body or the tail. 

In reference 9^ nxmnerical calculations were carried out to determine 
the positions of the vortices. at various distances behind the iTings of 
several plane' wing-body combinations . In order to make use of those 
calculations, we shall consider a specific slender ■t/ing-body combination 
and a specific tail length. Thus let 

— = 0.6 and ^ = 10 

Bq &o 

Also, to insure the validity of the slenderness assumption, we choose the 
case treated in reference 9 for A = 2/3, M = 2 where A is the wing 
aspect ratio and M is the free-stream Mach number. The body is cylindri- 
cal behind the wing trailing edge so that ax = Sq, and it will be assiamed 
that the vortex positions are not influenced by the tail. Thus from 
figure 10 of reference 9 


r = 1.27 Uoaso 


and the position of the vortex at the base plane x = I is given by 

= yi + izi where yi and Zx are obtained from figure 7(a)' of refer- 
ence 9 and Zx is transferred to body axes by adding cxd to the values 
in the figine. 

The above information can be put directly into equation (15) to 
determine the interference lift as a function of the ratio of tall span 
to wing span Sx/sq for all angles of attack for idiich reference 9 has 
supplied the vortex positions. The resiiltlng curves for the interference 
lift coefficient, based on the gross -tfing area, are shown in figure 1.^ 

It can be seen that the interference lift is nonlinear with respect to 
the angle of attack and becomes more negative as the ratio of tall span 
to wing span Increases. At sx/ Sq = 8 ./bq = 0.6, there is no tail and 
the Interference lift is the lift on the cylindrical afterbody. 


CONCLUDING REMARKS 

The analysis reported in reference 8 has been extended to permit the 
calculation of the total forces and moments (except drag) for slender 
^Thls calculation requires taking the real part of the square root 
indicated in equation (l5)- Since this process is not entirely straight- 
forward if ambiguities are to be avoided, the procedure is given in 
Appendix B. 
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•^•ring-'body-tail combinations of general cross section performing quasi- 
stationary maneuvers in a compressible fluid. Expressions have been 
developed in tbe present paper for the forces and moments due to the 
influence of ^rfng wakes’ on tbe other airplane components, and these eq.ua- 
tions were used to calculate the Interference lift of a plane 'vring-body- 
tail coniblnation. 

The essential q.uantlty required for the calculations is shown to 
be proportional to the impulse of each shed vortex and its image vortex 
in a transformed circle plane, and it is demonstrated that the calciilatlon 
of the interference lift in steady straight fli^t req.uires the determina- 
tion of this quantity only at the traiJLing edge of the wing and at the 
base of the configuration. 


Ames Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Moffett Field, Calif., Sept. 9, 1955 



APEEItDIX A 


LIST OF IMPORTANT SYMBOLS 

coefficient of term in Laurent eacpansion of the complex 

potential F(0 with no vortices in the field 
coefficient of term in expansion of the additional complex 

potential F’(0 due to shed vortices and their images 
body radius 

Uq dS 

coefficient of In t in expansion of F(^);B = 7 t — 3- 

ot ux 

distance from airplane nose to pivot point 

a function of x and t -vjhich contributes nothing to the forces 
and moments considered herein 

tall length, I - Xg^ 

con5)lex potential 9 + ii|f trith no vortices in the field 
additional complex potential due to shed vortices and their images 
length of airplane 

force in the z direction (approximately lift) 
rolling moment about the x axis 
pitching moment about pivot point x = ci 
number of external (free) vortices 
yaT'/lng moment about pivot point x = Ci 
angular rolling velocity about the x axis 
angular pitching velocity about the y axis 
V + iW 

angular yai'/lng velocity about the z axis 
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radiiis of transformed circle corresponding to airplane cross 
section 

cross-sectional area of the airplane 
local semispan of nlng or tail 
maximum semi span of wing (at x = xro) 
maximum semispan of tail (at x = Z) 
time 

component of flight velocity along negative x axis 

component of flight velocity along positive y axis 
(Vo = UoP if P = 0 ) 

Vo - r(x - Cl) 

component of flight velocity along positive z axis 

(Wo = -Uoo, if p = 0) 

¥0 - g(x - Cl) 

force in the y direction (side force) 

Cartesian coordinates fixed in the body (x rearward, y to 
starboard, z upward, origin at the body nose) 

distance from airplane nose to -vTing trailing edge 

angle of attack (angle het\-reen arbitrarily chosen xy plane 
and flight direction) 

angle of sideslip (angle bebi-reen xz plane and flight direction) 

circulation strength of kth external (free) vortex, positive 
counter clocto-Tis e 

fluid mass density 

y + iz 

complex coordinate of centroid of cross-sectional area 
position of kth external (free) vortex 
complex coordinate in transformed circle plane 
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a]j position of kth exteraal vortex in transformed circle plane 

position of kth external vortex relative to its imag e in the 

transformed circle plane, c% - 

Ok 

9 velocity potential 

Tjr stream function 


Special Notations 


/ 


subscript indicating vortex interference 

contour integral taken once round the body cross section in the 
positive (counterclockwise) sense 


X imaginary part 


R real part 

( ) complex conjiigate of ( ) 
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APFEHDIX B 


EVAIBTAIEIOB OB’ THE REAL PART OP A COMRLEX SQUARE ROOT 


A coniplex square root of the type - s^ ■where ^ = y + iz and 

s is real has hranch points at ^ = ±s and can he -written in its 
factored form 





- s)(£ + s) 

Roir each factor can he -I'/ritten in 
polar coordinates referred to one of 
the hranch points. Thus (see sketch 
(e)) let 

^ - s = rxs^^; ^ + s = rge^^ 

where 7 and A are hoth limited to 
a range of, say, 0 to 2jt (to give 
the proper sign changes throu^ the 
line segment shown) . This is commonly 
called a bipolar coordinate system 
and enables one to -write 


= *7rar2 



= Vrjra 


cos 


7 + A' 


+ i sin 




so that -the real part is 


B 


= JrxTz cos 


1-jhere 


Tz = 


= l*v/(y - s)^ + 
7(y + s)^ + 


7 = tan* 


-xf z 


y - s 


A = tan' 
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The sign of the real part of the square root is therefore determined hy 
the quadrants of 7 and A ■^^hich depend on the position of the point ^ 
relative to the branch points ^ = ±s and also on the range of 7 and 
A specified above . 
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